Half-filled Hubbard Model on a Bethe lattice with next-nearest neighbor hopping 
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We study the interplay between Neel-antiferromagnetism and the paramagnetic metal-insulator- 
transition (PMIT) on a Bethe lattice with nearest and next-nearest neighbor hopping t\ and t2. 
We concentrate in this paper on the situation at half-filling. For t2/t\ 1 the PMIT outgrows the 
antiferromagnetic phase and shows a scenario similar to V2O3. In this parameter regime we also 
observe a novel magnetic phase. 



I. INTRODUCTION 

Understanding correlation effects is one major goal of 
condensed matter physics. Strong correlations manifest 
themselves in various forms. The paramagnetic Mott- 
Hubbard metal-insulator-transition (PMIT)i is a well- 
known and interesting example. With increasing inter- 
action strength the Fermi liquid state breaks down at a 
critical value and an insulator is formed. 

Another fundamental example is magnetism, where 
electrons reduce the energetic cost of the Coulomb- 
interaction by ordering. Both effects can of course occur 
simultaneously and are the heart of the extremely rich 
phase diagram of e.g. transition metal compounds like 
for example V2O3 or LaTiOg^i^. 

Besides strong correlations, another major ingredient 
for the understanding of the phase diagram of compounds 
like V2O3 is frustration. V2O3 crystallizes in the corun- 
dum structure with the V-ions located on a honeycomb 
lattice in the ab-plane, while along the c-axis a more com- 
plicated coordination is observed, which induces frustra- 
tion of the magnetic interactionsi. Nevertheless does the 
phase diagram of V2O3 show an antiferromagnetic phase 
at temperatures below T/v ~ 180K. Upon doping with 
Ti one may suppress this order. Such a doping with a 
smaller ion can be viewed as internal pressurcii hence 
the suppression of the magnetic order is commonly in- 
terpreted as happening through an increase of the band- 
width respectively a decrease of the correlation effects. 
Consequently, the critical Ti doping is conventionally re- 
lated to the existence of a lower critical value of the elec- 
tronic interaction parameter. At higher temperatures 
the antiferromagnetic state becomes unstable towards 
a paramagnet and one can eventually observe a para- 
magnetic metal-insulator-transition up to temperatures 
T « 400 A'. 

Frustration is a quite common feature in real ma- 
terials. Very interesting examples for frustrated sys- 
tems are layered organic compounds like k-(BEDT- 
TTF)2Xi^>&i^^aodia2i3aia5a6i7^ They have a similar 

phase diagram as the high-temperature superconductors 
(HTSC)^. The phases of these organic systems are con- 
trolled by pressure and frustration rather than by dop- 
ing as in HTSC— . They are usually described by an 
anisotropic triangular lattice, and changing the anion (X) 
in these systems modifies the frustration of the lattice. 



Besides superconductivity also magnetic ordering and a 
PMIT can be found. 

These two examples by no means exhaust the zoo of 
materials showing such interplay or competition between 
PMIT and ordered phases^. For example, rare-earth 
compounds like Ce(Rh,Ir)i_2;(Co,Ir)2;In5 do show a simi- 
larly bizarre phase diagram^S. Besides their usually com- 
plicated lattice structure another challenge for a theoret- 
ical description of such compounds is that the presence 
of elements with partially filled d- or f-shells in principle 
requires a multi-orbital description to account for effects 
like Hund's or spin-orbit coupling properly. Furthermore 
the residual degeneracies in the solid state crystalline en- 
vironment lead to degenerate multiplcts which in turn 
can give rise to even more complex structures like orbital 
order or polaron formation (see e.g. Imada et al.— for an 
overview and references). 

Although all these ingredients play an important role 
for a quantitative theoretical description of transition- 
metal or rare-earth compounds, we here want to focus on 
the one-orbital situation, in particular on the relation be- 
tween PMIT and antiferromagnetism. This restriction to 
a simpler but by no means trivial situation will enable us 
to investigate the relation between these two paradigms 
of correlation effects with a small and controllable set 
of parameters and thus obtain some hint to how both 
phases interact. A model suitable for analyzing this kind 
of physics is provided through the Hubbard mode l^^'^^i^^ 
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where c\^{cia) creates (annihilates) an electron with spin 
a at site i and rtij (rti| ) is the density operator for spin 
up(down) at site i. The parameters represent the 
hopping amplitude from i to j and U is the interaction 
strength. In this paper we will measure the interaction 
relative to the bandwidth, which is related to the hop- 
ping amplitude. Although at first sight very simplistic, 
this model is highly nontrivial. Besides other methods, 
especially in one dimension, progress in understanding its 
physics was achieved by the development of the dynam- 
ical mean field theory (DMFT)^. The DMFT is a very 
powerful tool for analyzing strongly correlated lattice sys- 
tems, mapping the lattice problem onto a quantum impu- 
rity problem, which has to be solved self consistently. For 
solving this impurity problem for arbitrary interaction 
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strengths and temperatures we here use Wilson's numeri- 
cal renormalization group ^^i^^ . An interesting fact is that 
the only information about the lattice structure, which 
enters a DMFT self consistency calculation, is the local 
density of states (DOS) of the non-interacting system. 
We performed our calculations for a Bethe lattice with 
nearest-neighbor (nn) and next-nearest-neighbor (nnn) 
hopping ti and ^2, respectively. The DOS in this case 
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FIG. 1: (color online) DOS for increasing nnn-hopping t2. 
Observe the van-Hove singularity at the lower band edge. 
p{u}) and ij are scaled with the bandwidth W = 4t2 + 2ti -|- 
tl/iUi) [t2 > l/4ti). 

can be calculated using a topological ansata^li^. Start- 
ing from a particle-hole symmetric DOS at ^2 = the 
density of states becomes now asymmetric with increas- 
ing t2 (see Fig. [T]) and develops a van-Hove singularity 
at the lower band edge for positive and increasing nnn- 
hopping t2- In contrast to ^2 = 0, where the particle-hole 
symmetry can be employed to fix the filling at (n) = 1 
precisely, the asymmetry present for <2 7^ makes it more 
difficult to perform calculations with the filling kept at 
(n) = 1 with sufficient accuracy. Thus DMFT calcula- 
tions typically take much longer here due to the neces- 
sary adjustment of the chemical potential. Of course the 
Bethe lattice does not represent a lattice realized in real 
materials. However, in contrast to the hypercubic lattice 
with infinite coordination number the Bethe lattice has 
a compact support and thus possesses band edges, which 
provides a more realistic scenario. 

Since the early days of DMFT, there have been many 
contributions by different groups to the subject of the 
PMIT and antiferromagnet ism , ^"^i jjowevcr, frustra- 
tion effects up to now where introduced in DMFT typi- 
cally within the so-called two-sublattice fully frustrated 
mode l^"^i'^°'^^i'^^'^'^i'^^ , which results in a particle-hole sym- 
metric DOS even with frustration. As side effect, this 
way of introducing frustration leaves the paramagnetic 
phase unchanged. For the non-frustrated system the 
PMIT is then completely covered by the antifcrromag- 
netic phase, which exists for half-filling for all finite values 



of [7—1^. For the frustrated system, on the other hand, 
there exists a lower critical value for the interaction U, 
which increases with increasing frustration. It was fur- 
thermore found that the Neel-temperature decreases with 
increasing frustration such that the PMIT outgrows the 
antiferromagnetic phase^i. In early calculations using 
this way of introducing frustration based on exact diag- 
onalization studies of the two-sublattice fully frustrated 
model^ii^^iSj the authors also found parameter regions 
in the phase diagram where an antiferromagnetic metal 
appeared to be stable. However, this antiferromagnetic 
metal phase was later traced back to numerical subtleties 
in the exact diagonalization procedure and shown to be 
actually absent from the phase diagram^!. 

The first attempt to study the Hubbard model on 
the Bethe lattice with correct inclusion of nn- and nnn- 
hopping has been performed rather recentlj^. In this 
work the authors concentrated on the paramagnetic 
PMIT and found phase-separation between the insulat- 
ing and metallic phase. 

In this paper we investigate the PMIT as well as anti- 
ferromagnetism and concentrate on the competition be- 
tween the paramagnetic phase including the PMIT and 
the antiferromagnetic phase at intermediate and high 
grades of frustration. We especially look at the case 
t2 —>■ ti and raise the question, if the scenario of the out- 
growing PMIT, proposed in Zitzler et alM, stiU holds 
for the correct asymmetric density of states. The paper 
is arranged as following. After this introduction we start 
with a brief look at the PMIT, followed by a discussion of 
the phase diagram at half-filling including antiferromag- 
netism and the PMIT. The next paragraph addresses es- 
pecially the case of very strong frustration and the ques- 
tion how the magnetic order is realized there. The paper 
will be closed by a summary of our results and an out- 
look. 



II. METAL-INSULATOR-TRANSITION 

The metal insulator transition for the Bethe lattice 
with nnn-hopping has been analyzed by Eckstein ct al.— 
within the self-energy functional approach^. They par- 
ticularly focused on ^2/^1 = 3/7 and discussed an un- 
expected occurrence of phase separation in the param- 
agnetic state between a Mott-Hubbard insulator and a 
correlated metal at and near half-filling. Here we want 
to investigate the behavior of the system as function of 
increasing frustration. Due to symmetry there is no dif- 
ference between ^2 and — ^2- The calculations were done 
using Wilson's NRG as impurity solver for the DMFT. 
with A = 2, 1800 states kept per NRG step and a loga- 
rithmic broadening b = 0.8 to obtain spectral functions. 
We want to note at this point that the choice of NRG 
numerical parameters does not infiuence the qualitative 
nature of the results. We observe, however, small de- 
pendencies on A and b, which tend to become more pro- 
nounced close to phase transitions and may result in sys- 
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tematic errors in numerical values for critical parameters 
of the order of < 5%^ 

Figure [2] shows the paramagnetic metal insulator tran- 
sition for various values tijtx. As energy-scale we choose 
the bandwidth 



W 



4ti for < |t2| < ii/4 
4|t2|+2ti+t2/(4|t2|) for |t2| ><i/4 



of the non-interacting system. Note that these results 
are obtained by artificially suppressing an antiferromag- 
netic instability. We will come back to this point later. 
The occupation was kept fixed at n = f ± 0.005 by ad- 
justing the chemical potential. Note that in contrast to 
the case with i2 = it is not possible to achieve n = \ 
here within numerical precision due to the asymmetric 
DOS (see Fig. [J). For increasing ^2/^1 ^ 1 the PMIT 
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FIG. 2: (color online) The transition lines for the PMIT for 
different frustrations as function of temperature and interac- 
tion strength. For each frustration the right line represents 
the transition from the metal to the insulator, while the left 
line represents the transition from the insulator to the metal. 
Symbols mark the calculated data points, the lines are fits 
meant as guide to the eye. 



III. ANTIFERROMAGNETISM AT FINITE ±2 

We now allow for antiferromagnetic ordering in our 
calculations. To this end we reformulate the DMFT for 
an AB lattice structuro^i^ to accommodate the Neel 
ordering and initialize the calculation with a small stag- 
gered field, which is turned off after one DMFT itera- 
tion. The system then either evolves into a paramag- 
netic or antiferromagnetic state with increasing number 
of DMFT iterations. Figure [3] shows the resulting phase 
diagrams for t\ = 0.6 (upper panel) and tijtx = 0.8 
(lower panel) for different temperatures and interaction 
strengths. The small black points show the locations. 
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is shifted towards lower interaction strengths and also 
lower temperatures. While the shift in the interaction 
strength is rather moderate, we notice a large difference 
in the temperature of the critical endpoint between the 
unfrustrated and highly frustrated system. This observa- 
tion of course renews our interest in the question, to what 
extent long-range hopping can help to push the param- 
agnetic MIT out of the expected antiferromagnetic phase 
for reasonable magnitudes of to create a phase diagram 
similar to the one found for V2O3. The scenario proposed 
by Zitzler et al.— relied on the fact that the paramagnetic 
phase largely remains unaltered with increasing ti- As 
the Neel-temperature for the antiferromagnet is reduced 
at the same time, the PMIT can eventually outgrow the 
antiferromagnetic phase. 



FIG. 3: (color online) The upper (lower) panel shows the 
T — U phase diagram for t2/ti = 0.6 (0.8). The shaded area 
represents the antiferromagnetic phase, while the white area 
represents the paramagnetic phase. The lines show, where the 
PMIT in the paramagnetic phase would occur. The points 
denote the parameter values, where DMFT calculations were 
done. From these points the shaded area was constructed as 
guide to the eye. Additional calculations were performed to 
find the PMIT-lines. 

where calculations have actually been performed. From 
these data the shaded areas were constructed represent- 
ing the antiferromagnetic phases. This of course means 
that the phase boundaries shown here must be consid- 
ered as guess only. However, as we do not expect any 
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strange structures to appear, this guess will presumably 
represent the true phase boundary within a few percent. 

The full lines in Fig. [3] are the PMIT transitions. Note 
that for both diagrams the same division of axes was 
chosen. 

In contrast to the Hubbard model on a bipartite lattice 
with t2 = 0, there now exists a finite critical value U^^ ^ 
below which no antiferromagnetism can be stabilized 
even for temperature T 0. With increasing frustration 
the paramagnetic-antiferromagnetic transition is shifted 
towards higher interaction strengths and lower tempera- 
tures, while the PMIT is shifted towards lower interac- 
tions strengths. So obviously the PMIT is shifted towards 
the phase boundaries of the antiferromagnetic dome. So 
far this is the expected effect of the nnn-hopping which 
introduces frustration to the antiferromagnetic exchange. 
However, note that although ^2/^1 = 0.8 represents al- 
ready a very strongly frustrated system, the PMIT still 
lies well covered within the antiferromagnetic phase. 

Let us now have a closer look at the paramagnetic- 
antiferromagnetic transition. Here, Zitzler et al.— made 
the prediction that one has to expect a first order transi- 
tion close to the critical U^^ at low temperatures; while 
at larger values of U again a second order transition was 
found. 
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transitions from the antiferromagnetic to the paramag- 
netic state with decreasing interaction strength. In the 
upper panel (small U) one can clearly see a hysteresis of 
the antiferromagnetic transition. This hysteresis as well 
as the jump in the magnetization are clear signs for a 
first order transition. This antiferromagnetic hysteresis 
is very pronounced for strong frustration but numerically 
not resolvable for example for 12/ ti = 0.2. We believe 
that the hysteresis region shrinks with decreasing t2 and 
eventually cannot be resolved anymore with numerical 
techniques. The whole temperature depending hysteresis 
region can be seen in Fig. [5]for the case ^2/^1 = 0.8. One 
can see clearly the shrinking of the hysteresis region with 
increasing temperature. Note that such a hysteresis is 
also found in the two-sublatticc fully frustrated model^ 
which means that this quite likely is a generic effect in 
frustrated systems at intermediate coupling strengths. 

The lower panel in Fig. 0] shows the staggered mag- 
netization for temperatures just below the correspond- 
ing Neel-temperatures and at large interaction strengths. 
Here the magnetization vanishes smoothly, which is the 
behavior expected for a second order phase transition. In 
summary we thus find a first order transition at the crit- 
ical interaction U^^ where antiferromagnetism sets in, 
and a second order transition for the large Coulomb pa- 
rameter U ^ W. The merging from both transition lines 
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FIG. 4: (color online) Staggered magnetization versus interac- 
tion U/W for two different temperatures and 12/ ti = 0.8. In 
the upper panel there are for each temperature two transition 
lines, representing either increasing or decreasing interaction 
strength. The region between both lines embodies a hysteresis 
region. The lower panel shows the transition for large inter- 
action strengths. Here no hysteresis region could be found, 
but a smooth transition. 

Figure |4] shows the staggered magnetization for dif- 
ferent temperatures and interaction strengths at fixed 
^2/^1 = 0.8. The upper panel collects data for the tran- 
sition at low temperatures at the lower edge of the an- 
tiferromagnetic phase. The full lines represent the tran- 
sition from the paramagnetic to the antiferromagnetic 
state with increasing interaction strength for two differ- 
ent temperatures, while the dashed lines represent the 



FIG. 5: (color online) Phase diagram for t2/ti = 0.8 including 
the temperature depending hysteresis region (between dashed 
and full line). The meaning of symbols is as in Fig. [3] The 
shaded area is as before meant as guide to the eye. 

is an interesting point in itself. There must be a critical 
point where the first order transition changes into a sec- 
ond order transition. It is however not possible to resolve 
this merging within DMFT/NRG. First, the logarithmic 
discretization of the temperatures within the NRG does 
not allow to resolve this merging-region with arbitrary 
precision. Second, the magnetization of the system be- 
comes very small in this region, so it is not possible to 
distinguish between a (tiny) jump and numerical artifacts 
of a smoothly vanishing order parameter. Consequently, 
we cannot decide anymore of what order the transition 
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will be. 

Antiferromagnetic metallic phases at half-filling were 
reported in earlier publications^i^i^^i^. In our calcula- 
tions we saw no evidence for an antiferromagnetic metal- 
lic state at half-filling. Especially for strong frustration 
12/ ti K, 0.8 the system directly jumps from a param- 
agnetic metallic solution into an antiferromagnetic insu- 
lating solution with high magnetization. In the papers 
cited, the region showing an antiferromagnetic metallic 
solution broadens with increasing ^2- This prediction we 
clearly cannot confirm, as discussed above. Only in sys- 
tems with small to intermediate frustration there are nar- 
row interaction regimes where we observe a small finite 
weight at the Fermi level. One must however consider 
that the occupation number is not exactly one but only 
within 0.5%. Also it was sometimes difficult to stabilize a 
DMFT solution in these regions. In summary, we cannot 
see any clear signs for an antiferromagnetic metallic state 
at half-filling in our calculations. If any exists, then only 
for rather low frustration in a very small regime about 
the critical interaction. To what extent these rather spe- 
cial conditions can then be considered as realistic for real 
materials is yet another question. 

IV. NEARLY FULLY FRUSTRATED SYSTEM 

In this last paragraph we want to study the situation, 
in which ti and ^2 arc comparable in strength. Interest- 
ingly, there has been no attempt to calculate the phase 
diagram on a mean-field level in the strongly frustrated 
model t2 ~ ti. Therefore, before discussing the results of 
the DMFT calculations for strongly frustrated systems 
h/ti ~ 1 let us try to gain some insight into the physics 
we must expect, by inspecting classical spins on a Bethe 
lattice with nn- interaction Ji and nnn- interaction J2- Al- 
lowing that nearest-neighbor spins enclose an angle 6 one 
ends up with the energy functional 

E/2N = JiZcos{e) 
z-i 

+ J2Z (cos(6i)2 + sin(6')2 cos(27rV^)) (2) 

i=l 

Performing the same limits and scaling as in DMFT one 
finds (see Appendix A) 

E/2N ^ Ji cos(6') + J2 cos{ef. (3) 

Thus, the Neel-state with 6* = tt is the stable ground 
state for J2 < ^Ji, while one finds a spin wave with 
6* = TT — arccos( J1/2J2) for J2 > jJi- 

For the DMFT calculations we can allow only for solu- 
tions commensurate with the lattice. This however will 
possibly be inconsistent with the spin structure favored 
by the system. If, for example, we perform a calculation 
focusing on the ferromagnetic solution within a param- 
eter regime, where the system wants to order antiferro- 
magnetically, DMFT will not converge. To investigate 



spin wave states with periodicities with more than two 
lattice sites, one has to set up the correct DMFT self- 
consistency equations respecting the lattice structure. 
While for a system on an infinite Bethe lattice with nn 
hopping only it is straightforward to extend the DMFT 
to commensurate magnetic structures with periodicities 
of more than two lattice sites, we did not succeed in devis- 
ing a scheme that allows for such calculations for systems 
with nnn hopping. The reason is that one has to parti- 
tion the lattice into an ABCD . . . structure. However, 
the nnn hopping makes it impossible to uniquely identify 
the connectivity of the respective sublattices. A method 
proposed by Fleck et al.— for the two-dimensional cubic 
lattice is not applicable in our case. 

We thus only allowed for paramagnetic, ferromagnetic 
and antiferromagnetic solutions in our calculations. The 
resulting phase diagrams for t2 ^ ti are shown in Figs, 
ini and [71 Figure [H] displays the ground states for differ- 
ent grades of frustration and interaction strengths. For 
ti/ti < 0.95 the phase diagram has the same structure 
as for small and intermediate ^2- The critical interac- 
tion strength U^^ necessary to stabilize the Neel-state 
increases and for all values above U^^ we find an anti- 
ferromagnetic phase with a hysteresis region at the phase 
boundary. For 0.95 < 12/ ti < 1 the critical value U^^ 




FIG. 6: (color online) Ground state (T = 0) phase diagram 
for different strengths of frustration as function of the interac- 
tion. The brindled regions are hysteresis regions for increas- 
ing or decreasing interaction. The phase boundaries of the 
incommensurate (IC) phase are only qualitative (explanation 
see text). 

one needs to stabilize the Neel-state increases dramati- 
cally. For ^2 = ti finally we do not find an antiferro- 
magnetic Neel-state at all for any interaction strength 
U. Our DMFT calculations however indicate that in this 
range of 12/ ti there actually does exists another mag- 
netic phase. Namely, for sufficiently small temperatures 
one obtains a finite spin polarization in every DMFT it- 
eration. However, the DMFT does not converge to a 
unique state as function of DMFT iterations (see also 
Fig. [HI). In the phase diagrams in Figs. [5] and [7| we 
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U/W 



FIG. 7: (color online) Phase diagram T - f/ for ta/ii = 0.98 
including the PMIT. The brindled area represents again the 
hysteresis region. 




DMFT iteration 

FIG. 8: Example of a non-convergent DMFT calculation. The 
figure shows the staggered polarization over the DMFT iter- 
ation number for t2/ti = 0.98, U/W = 1.35 and T « 0. The 
lines are meant as guide to the eye. 

have named this regime the incommensurate phase (IC 
phase) . In this parameter regime the Neel-state becomes 
unstable towards the behaviour shown in Fig. [S] Here 
one can switch between a conventional Neel-state and 
the IC phase by only a small change of the interaction 
strength. Note, that the pliase boundaries sliown in the 
figure must be taken with some care as we cannot com- 
pare the energies of the Neel-state and this IC phase to 
properly determine the phase boundaries. As we observe 
precisely the same behavior for all investigated values 
t2/ti {0.96, 0.97, 0.98, 0.99, 1.0} we are convinced that 
the ground state in this region is an incommensurate 
state, as to be expected from our results for S — oo. 
Similar observations also hold for finite temperatures as 
shown in Fig. [71 where the T — U phase diagram for fixed 
^2/^1 = 0.98 is displayed. For increasing interactions and 
T = there first is a transition from a paramagnetic 
metal to the IC phase and for U/W w 1.6 from the IC 



phase to the Neel-state. For increasing temperature the 
IC phase eventually becomes unstable towards the Neel- 
state. In Fig. [7] one can also see the PMIT lines. As one 
can see it lies within the hysteresis region of the magnetic 
phases but clearly outgrows both magnetic phases. This 
is the scenario described in Zitzler et al.— . 



V. SUMMARY 

We studied the DMFT phase diagram of the Hubbard 
model at half-filling in the presence of nn- and nnn- 
hopping. In contrast to previous investigations we did 
our calculations for a Bethe lattice with proper nnn hop- 
ping t2, introducing a highly asymmetric DOS already 
for the non-interacting system. 

The first important observation concerns the paramag- 
netic metal-insulator transition, which is suppressed by 
increasing t2 , but at the same time shifted to lower values 
of the Coulomb interaction. 

The at t2 = ubiquitous antiferromagnetic phase 
on the other hand is suppressed up to a critical value 
U^^{t2) with increasing t2, as expected. Furthermore, 
a hysteresis region between the paramagnetic metal at 
small U and the antiferromagnetic insulator at large U 
develops, showing that the transition is of first order. 
Note that we did not observe any evidence for an antifer- 
romagnetic metal close to the phase boundary, nor did 
the PMIT reach out of the antiferromagnetic insulator 
up to values t2 = O.Sti. 

Thus far the observations are similar to the results 
found by Zitzler et al.— for the two-sublattice fully frus- 
trated Bethe lattice^!. The shift of the PMIT to lower 
values of U together with a moderate suppression of the 
critical temperature for larger t2 however motivated a 
more detailed investigation of the region of larger t2 ■ A 
simple argument based on classical spins with competing 
interactions showed that one has to expect an additional 
incommensurate phase here. In fact, as already anti- 
cipated qualitatively by Zitzler et al.— , for frustrations 
0.96 < t2/ti < 1.0 we eventually found that the PMIT 
lies within the hysteresis region of the antiferromagnetic 
phase for T = 0, but outgrows it in temperature. For 
such strong frustration we also found evidence for another 
magnetic phase besides ferromagnetism or antiferromag- 
netism. Unfortunately this phase could not be stabilized 
within our DMFT calculations, so its real nature remains 
open. In connection with our argument based on classi- 
cal spins we believe that we can interpret the observed 
structure as an incommensurate phase. This conjecture 
is further supported by the fact, that for ti = t2 we found 
no antiferromagnetic solution of the Neel-type, but only 
this frustrated magnetic phase. 

Especially the latter findings make it highly desirable 
to set up a scheme that allows to study commensurable 
structures with period beyond Neel-type for arbitrary lat- 
tice structures including longer-ranged hopping. 
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APPENDIX A: CALCULATION FOR 
VECTOR-SPINS 

Here we want to present the calculation for 3- 
dimensional vector spins on a Bethe lattice with antifer- 
romagnetic coupling between nn- and nnn-lattice sites. 
We take Z nearest neighbors and interaction strengths 




FIG. 9: (color online) Left: Bethe lattice Z = ?> with sites 
numbered according to the vector spins on the right. The 
nearest neighbors of site 1 must lie on a circle so there is an 
angle d between 1 and 2. Similarly the nearest neighbors of 
one of the 2-spins must lie on a circle including the 3-spins 
and the 1-spin. 



Ji between nn-sites and Ji between nnn-sites. The last 
parameter entering this calculation is the angle Q be- 
tween nn-spins. Although the initial assumption that two 
neighboring spins form an angle theta may seem some- 
what restrictive, we are not aware of other configurations 
with lower total energy^. We want to minimize the en- 
ergy with respect to this angle. According to Fig. [9] 
the nn-spins of one spin, must lie on a circle. The spins 
ending on the circle are all nnn-spins. Due to the anti- 
ferromagnetic interaction J2, we assume that they want 
to maximize the angle between them. Since there are Z 
spins on each circle, we assume they will have angle I-k jZ 
projected on the circle. Using now simple trigonometry, 
the angle between nnn-spins is given by 

2-(2i?2_2i?2 cos(27rV^)) 

cos 7 = 

' 2 

cos(6l)2 4-sin(6l)2cos(27ri/Z), 

where i runs from to Z — 1, giving the different positions 

on one circle. Inserting this into the Hamiltonian 

E = Ji SiSj + J2 SiSj, 

one finds for the energy 
E/2N = JiZcos{e) 

+ J2Z {cos{ef + sin(6l)2 cos{27Ti/Z)) . 

Performing now the limit Z — > cx) and scaling JiZ 

and J2Z J2/Z one finally obtains for the energy per 

lattice site 

Ez=oo{0)/{'2-N) = J* cos{e) + J* cos{0f . 

One can now see that the Neel-statc 6* = tt is stable for 
J2/J1 < 1/2, because (fE{e = TT)/de^ = J* - 2J*. 
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